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Abstract. It is difficult to derive the solid-fluid transition from microscopic models. We introduce 
particle systems whose potentials do not decay with distance and calculate their partition function 
exactly using a method similar to that for lattice systems with infinite-range interaction. In particular, 
we investigate the behaviors of examples among these models, which become a triangular, body-centered 
cubic, face-centered cubic, or simple cubic lattice in low-temperature phase. The transitions of the first 
three examples are of the first order, and that of the last example is of the second order. Note that 
we define the solid phase as that whose order parameter, or Fourier component of the density, becomes 
nonzero, and the models we considered obey the ideal-gas law even in the solid phase. 


1. Introduction 

A fundamental problem in condensed matter physics and statistical physics is the theoretical derivation 
of solid-fluid phase transitions from microscopic Hamiltonians. There are several ways to define the 
solid and fluid phase. Some researches distinguished two phases from each other by the existence of the 
shear modulus [T]. Other researches considered the solid phase as the state with non-uniform periodic 
density. This definition is originated back to the theories by Vlasov[^, and also Kirkwood-Monroej^. 
In the standpoint of modern statistical physics, the latter definition is interpreted as the continuous 
translational symmetry breaking and the emergence of the corresponding order parameters. In this 
paper, we adopt the latter in order to focus on the symmetry breaking and the order parameters. 
Although many efforts have been made to develop the theory of solid-fluid transitions, most of them 
are based on phenomenological theory, and some do not take the translational symmetry breaking into 
account explicitly[31 01 [S]. Only a few theoretical approaches, including density functional theory, have 
successfully overcome these problems to some extent [iiili. Hence, to more deeply understand solid- 
fluid transitions, it is meaningful to search for and consider exactly solvable models that exhibit solid-fluid 
transitions. One of a few examples in which the solid-fluid transition is derived exactly is the research 
of Carmesin and Fan m- They considered a one-dimensional model with a cosine potential and hard 
core potential, and used the Hubbard-Stratonovich transformation for the cosine potential, as in the 
infinite-range XY model. However, their research used a special method to treat the hard core potential 
and was difhcult to apply to more general cases. Therefore, it is natural for us to try to simplify their 
model by removing the hard core potential and to extend the discussion to higher-dimensional cases. 

In this paper, we propose higher-dimensional models whose particles interact with each other only by 
the summation of the cosine potentials and derive the solid-fluid transition exactly. First, we introduce 
the models in section 2, and then we calculate the partition function in sections 3 and 4. We next 
investigate the behavior of the order parameter and the order of the phase transition in section 5, and 
finally summarize this paper in section 6. 

2. Models 

We consider a d-dimensional model depending on the values of fe^, or the smallest reciprocal lattice 
vectors of the crystal under our consideration: 

i ol i,j a i 
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where N is the number of particles; Xi and are the position and momentum of the tth particle, 
respectively. We treat the spontaneous symmetry breaking by using the method of Bogoliubov quasi¬ 
averages, namely weak external field ha with positive values is introduced to break symmetry. The set 
of ka values contains all the reciprocal lattice vectors with the smallest sizes. When we introduce the 
ka values, we regard a pair of reciprocal lattice vectors that are —1 times each other as one vector. For 
example, the set of ka values is given as 
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for the d-dimensional simple cubic lattice. Here, k is the length of ka- Notice that we introduce 
Hamiltonians which are different from each other for each crystals we consider. We consider crystals 
which include only one position for particles to rest in their primitive cells and span their reciprocal 
lattice space by ka- Hence, we do not consider crystals such as hexagonal lattices and hexagonal close- 
packed lattices, which have several positions for particles to rest. 


3. Calculation of the partition function 

In this section, we calculate the partition function of these models: 

' i '' i 

' fij 


■ exp 
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EE cos ka ■ {Xi — Xj) -\- fi ^ ha ^ COS ka ’ Xi 
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First, the cosine potential is transformed into one-body potentials with auxiliary variables qca, Qsa using 
the Hubbard-Stratonovich transformation: 


exp 
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Then, the partition function is expressed as 
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This method of calculating the partition function is similar to that for the infinite-range XY model and 
that used in previous works for one-dimensional systems lano]. We next maximize the right-hand side 
of equation ([9]) under the condition that the values of = \/qfa~\-q^ are fixed. The maximum occurs 
when qca = ra, qsa = 0 because of the weak external fields {ha}. The validity of this analysis is discussed 
in the Appendix. Considering this discussion, the partition function is expressed as the maximum of an 
r„ function: 
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Taking the periodicity of the potential into account, we express the integral in equation (1101) as 

/ dx • exp < V] {2j5Jra cos feq • xj i = ({rq} , {fea}), 

Jn \ a ] Vp 

where Vp and $p represent the integrals over a primitive cell of the supposed crystal and are given as 
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The solid-fluid transition is investigated using this equation. At the end of this section, we discuss 
the thermodynamic properties of this system. The order parameters, or the Fourier components of the 
density, are given as 



J _^ 

/3IV dha 


\ogZ 
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(15) 


where fa {{ha}) represents the value of that yields the maximum in equation (1141) . The equation of 
state is derived from the partition function as 


d /logZ\ _ N 

W V /3 ) ~W' 


(16) 


This is nothing but the ideal-gas law, which holds even in the solid phase. This strange property results 
from the fact that although the cosine potentials tend to restrict the positions of the particles to a set of 
lattice points, it does not matter on which of these points particles rest. 


4. Analytic expression of the partition function 

In this section, we discuss an analytical expression of $p starting from equation First, the 

exponential in the integral is expanded as follows: 
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We ignore the values of the weak external field ha in this section because their role ended when they 
gave the value of the order parameter under the condition that ha —>■ -(-0 in equation (|15|) . The integral 
on the right-hand side of equation m is calculated as 


dx exp < i ''^^{ua — n'a)ka ■ x \ = 
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Consequently, the contribution of the terms that have nonzero values of X)a(^c( ~ n'^)ka is sufficiently 
small. Hence, the partition function is given as 
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This is a general expression of the partition function, independent of the set of ka values. 

In the following, the partition function for each lattice structure is discussed. For the triangular 
lattice, taking the fact that fci -|- fc 2 + ^3 = 0 into account, we simplify equation (ITOl) as 
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where Ii, is the i^th modified Bessel function: 
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Assuming that this system does not break the symmetry under the permutation of a, and hence all r, 
have the same value r, we obtain 
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Although the case of the bcc lattice is more complicated, a similar expression for the partition function 
appears after some calculations: 
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For a simple cubic lattice, because the smallest reciprocal lattice vectors ki,k 2 , ■■■,kd are linearly 
independent, the nonzero terms in equation p9p fulfill Ua = rd^ for all a values. Hence, the partition 
function is written as 
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This is proportional to the dth power of the partition function of a one-dimensional system with the 
cosine potential. 


5. Behavior of the phase transition 

5.1. Temperature dependence of the order parameter 

In this section, we investigate the behavior of the order parameter r^ by calculating equation (HI 
numerically. If we examine equation (|19p in the previous section carefully, we find that this partition 
function is independent of k, or the length of ka. Hence, we let fc = I and also assume that all Tq have 
the same value as in the previous section. The results are shown in figures UHll To identify the order of 
the transition, we investigated the r dependence of the logarithm of exp (^—n'fiJr^) $p/I^, which appears 
on the right-hand side of equation (ITTl) . This function corresponds to —/3 times the Landau free energy 
except for the difference in the constant. For the triangular lattice, this function has two peaks, the 
heights of which are interchanged at the transition point, as we show in figure [5] Hence, the solid-fluid 
transition is a first-order transition in this case. The transition of bcc and fee lattices is also found to be 
a first-order transition by a similar investigation. However, the transition of the simple cubic lattice is 
of the second order because the behavior of the order parameter is the same as in the infinite-range XY 
model. Note that although the power series expansions of the partition function of the fee lattice model 
do not include odd-order terms of the order parameter, the transition in this model is of the first order. 
This is caused by the sixth- or higher-order terms. 
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Figure 1. Temperature dependence of the order parameter for the triangular lattice model 



Figure 2. Temperature dependence of the order parameter for the bcc lattice model 



Figure 3. Temperature dependence of the order parameter for the fee lattice model 
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Figure 4. Temperature dependence of the order parameter for the simple cubic lattice model 



Figure 5. r dependence of —/3F for the triangular lattice model at transition point without difference 
in constant 
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5.2. Effect of additional potential with higher wavenumber on the simple cubic lattice 


As we saw in section 15.11 the model of the simple cubic lattice has a second-order transition if {fee} 
includes only the smallest reciprocal lattice vectors regardless of the spatial dimension d. However, it 
is possible that the transition is of the first order if {/Cq} includes other reciprocal lattice vectors. We 
consider the case of d = 1 that contains k and 2k in {fea} as an example: 

2 j 

^{kc} ~ ^ 2^-^ {cosk{xi — Xj) -\- cos2k{xi — Xj)} — h''^coskxi. (27) 

i i,j i 

The partition function of this model is calculated by a method similar to that in sections |3] and |4] as 
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. (28) 


Consequently, the solid-fluid transition in this model is of the first order. The temperature dependence 
of the order parameter is shown in figure [Bl If we ignore the terms with a ^ 0 in equation (1281) , the 
summation appearing on the right-hand side of this equation reduces to the simple product of the zeroth 
modified Bessel function, and the transition is judged to be of the second order, like that of equation 
(IM)l . In other words, the terms with a ^ 0 cause the transition to be of the first order. 



Figure 6. Temperature dependence of the order parameter for equation isg 
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6. Summary 

We introduced two- or higher-dimensional models in which particles interact with each other by the 
summation of the cosine potentials and calculated their partition function. In particular, we investigated 
four examples that form a triangular, bcc, fee, or simple cubic lattice below the transition point. The 
phase transitions of the first three examples are the first order, and that of the simple cubic lattice is a 
second-order transition. 

Because the potentials of these models do not decay with distance and are treated similarly to 
the infinite-range model of spin systems, these models are thought to be called “mean-field models” 
of the solid-fluid transition. Hence, these models are expected to be used to constitute the “mean- 
field approximation” of solid-fluid transitions of simple particle systems, such as the three-dimensional 
Lennard-Jones system. In order to construct the approximation properly and obtain a non-trivial equation 
of state, not the ideal-gas law, it is necessary to introduce an effective coupling between the order 
parameter of the solid phase and the density or the volume of the system. 
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Appendix A. Proof of = ra,qsa = 0 in section [3] 


To prove that the integral in equation m reaches its maximum when gca = 'Ta, Qsa = 0 under the 
condition that qca + Qsa fixed, we first expand the integral in equation ([9|) in the same way 

as that in equation (HU. 
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By substituting (jA.ll) into (ID), the partition function is given as 
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The proposition we try to prove here will be fulfilled if the real part of each term appearing in the 
summation of equation (IA.3I) becomes maximum when qca = rg^qsa = 0. We show this fact starting 
from the following inequality; 
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-'■“ + 27 


(A.4) 


The equality is satisfied when q^a = ra,qsa = 0; hence, the integral in equation 0 takes its maximum 
in this case. Furthermore, considering the term when Ua' = n'^, = 1 for a = a' and na = n'^ = 0 for 
other a values, the inequality (IA.4I) becomes 


ReZc^'Za' = |7:|^ =rl+ 2qca • ^ 

< + Ij) . (A.5) 

Both sides of the inequality (IA.5I) are equal to each other if and only if qca = r^^qaa = 0. Because of this 
example, the case that qca = Ta , qsa = 0 turns out to be the unique one that yields the maximum value 
of the integral in equation dH). 
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